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Abstract 

Although conservative Hamiltonian systems with constraints can be formulated 
in terms of Dirac structures, a more general framework is necessary to cover also 
dissipative systems and gradient systems with constraints. We define Leibniz-Dirac 
(LD) structures which lead to a natural generalization of Dirac and Riemannian 
structures, for instance. From modeling point of view, LD structures make it easy to 
formulate implicit dissipative Hamiltonian systems. We give their exact characteri- 
zation in terms of bundle maps from the cotangent bundle into the tangent bundle. 
Their behavior under push-forward maps is also considered. Physical systems which 
can be formulated in terms of LD structures are discussed. 
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1 Introduction 

Dirac structures embody a number of geometric structures such as symplectic, Pois- 
son, foliation, complex geometries [TJ [2]. Since their first introduction [31 H] there 
have been a great number of work, which is still growing, have been done over the 
years. One of the most striking features of Dirac structures is that they can give a ge- 
ometric picture of Hamiltonian systems with constraints, holonomic or nonholonomic 
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[5j[6]. Nevertheless, Dirac structures are insufficient in formulating non-conservative 
Hamiltonian systems such as gradient systems or dissipative systems. In that sense, 
recently some attempts have been done to put these systems into a rather Hamilto- 
nian form. 

For example, in a generalization of Dirac structures is given in terms of an 
inner product of split sign on the Pontryagin bundle instead of the natural symmetric 
pairing. We specify this definition in order to cover the physical examples which are 
aimed to be put into the Hamiltonian context. In [3J, the authors use the notion of 
Leibniz structures [9] which is a generalization of Poisson structures, whose tensor is 
not necessarily skew-symmetric. Our approach is quite similar but we both work on 
the Pontryagin bundle and with systems with constraints on the manifold. In |10j . 
dissipative Hamiltonian systems with constraints are studied with Dirac's original 
method of reduced brackets. For other recent work on the generalizations of the 
conservative Hamiltonian systems we refer to [T2l [T3] and the references therein. 
Our motivation in this paper is to give a generalization of Dirac structures and study 
their geometric features in order to construct a general framework of non-conservative 
Hamiltonian systems. 

We define Leibniz-Dirac structures (LD structures in short) by weakening the 
properties of Dirac structures as follows: Let V be a vector space with its dual 
denoted by V * . A subspace V C V © V * is called a Dirac structure if it is maximally 
isotropic under the symmetric pairing 

({Vi,vl), (V2,V%))+ = (V*\V 2 ) + (V^l) 

for all (vi,Vi),(v2,v%) G V © V*, where (|) denotes the natural pairing between 
vectors and co- vectors. As a result, it is shown in [1] for a Dirac structure that the 
equations, which are called characteristic equations, 

p{V)° = VC\V* and p*(V) = {V n V)° 

are satisfied, where p and p* denote the projections from V © V* onto the first and 
second factor respectively, and (°) stands for the annihilation operator. Accordingly, 
there exist skew-symmetric linear maps uj : E = p(T>) — > E* and J : p*(T>) — > p*(T>)* . 
As the systems we aim to formalize can be given in terms of linear maps of the form 
V* — > V, we assume the second characteristic equation as the defining conditions of 
the generalizations of Dirac structures. They are called Leibniz-Dirac structures (LD 
structures for short). 

LD structures share some properties of Dirac structures such as being a La- 
grangian subspace for a suitable inner product. Their extensions on manifolds is de- 
fined the same way as in the Dirac structures. Several properties of linear and smooth 
LD structures are discussed. It is shown that LD structures are general enough to 
express possibly dissipative implicit Hamiltonian systems with constraints. From the 
discussion above, it is seen that LD structures can be pushed forward but cannot be 
pushed back, in the sense of |14j . We define LD maps which transform LD structures 
to LD structures. 

The paper is organized as follows. In Section [2] we define LD structures on mani- 
folds and give their characterization in terms of vector bundle maps from the cotan- 
gent bundle into the tangent bundle with Theorem [TJ Then it is shown that LD 
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structures are Lagrangian subspaces with respect to suitable symmetric pairings. 
Smooth LD structures on manifolds are defined in Section [TU where we also relate 
LD structures to Leibniz structures . In Section [3] we study admissible functions on 
manifolds with LD structures, then we study Hamiltonian dynamics of LD manifolds. 
We give present several physical examples which are also given in different formalisms. 
Maps that carry LD structures from one manifold to another manifold are studied 
in Section 14.21 The paper ends with some conclusions and future questions. 



2 Leibniz-Dirac structures 

In this section we develop the basic geometry of LD structures. First, we study linear 
LD structures and then extend it to smooth LD structures. Most of the results and 
notions are a refinement of the Dirac case pQ. 

2.1 Linear Leibniz-Dirac structures 

Let V be an n-dimensional vector space and V* be its dual space. Consider the 
direct product space V (B V* and denote the projections on V and V* by p and p* 
respectively. Also denote by W° the annihilator of a subspace W C V. 

Before proceeding in definition of linear LD structures, we give the following 
elementary but key result on subspaces of V V* pQ . 

Proposition 1. LetU be subspace ofV@V*, then the following equation holds: 

dim{U HV) + dim{p*{U)) = dim{U). (1) 

Proof. Let d\m{U) = k. If one chooses a basis for U, then this is equivalent to 
giving two linear maps a : R fc — > V and b : M. k — > V* such that the basis becomes 
(a ei , b e\ ) , . . . , (a e n , b e n ) . Since T> is fc-dimensional, 

kerankei-6 = {0} (2) 

is satisfied. Observe that T>DV = a(ker b) and p* (T>) = range b. But, dim(a(range b)) = 
dim(rangeb) by ([2]). Then by the fact that dim(a(ker b)) + dim(range&) = k it is 
seen that Equation ([U is justified. □ 

Definition 1. Let V be a vector space. A vector subspace T> C V © V* is called a 
Leibniz-Dirac structure (LD structure for short) on V if 

p*{v) = {vnv)° (3) 

is satisfied. 

Apparently, this is a generalization of Dirac structures. This and the other special 
examples will be clear after the following theorem. But we first give the following 
conclusion which plays an important role in what follows. 
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Proposition 2. LetT> be a a subspace ofV®V* with dim(V) = n, then the following 
are satisfied: 

(i) If D is a LD structure, then dim(T>) = n. 

(ii) If V is n- dimensional and p*(T>) C (Pfl V)° , then V is a LD structure on 

V. 

Proof, (i) follows from Proposition^ since dim((DnV r )°) = n — dim(T>nV). To show 
(ii), observe that p*(T>) C (Pfl V)° implies p*(V) = (Dfl V)° again by Proposition 

m □ 

Theorem 1. IfDisa LD structure on a vector space V , then there exists a linear 
map C : p*(V) ^V/(T>n V) whose kernel isVnV*. 

Conversely, given a subspace E C V and a linear map C : E° —> V/E, one can 
define a LD structure on V by 

V = {(v,v*);v = Cv*(modE)}cV®V*. (4) 

We first give the following lemma which is essential for the proof. It was originally 
given for Dirac structures in [TJ, and for LD structures the result was used in |7J 
without proof. 

Lemma 1. Let V be a n-dimensional vector space and V* denote its dual space. A 
pairing on V ®V* given by 

< («i,vj),(«2 > u5)»= (vt\v 2 ) + (v* 2 \v x ) -2m{v{,v* 2 ), (5) 

for all (v\, v I ), fajVri) G V @ V* where ^ : V* x V* — > R is a symmetric bilinear 
form, defines a non degenerate symmetric bilinear form of split sign. 

Proof. After choosing a proper basis for V, the result will be clear. First define the 
associated linear map ip : V* — > V to ^ by 

{v\\i)v* 2 ) =^(v\,v* 2 ) foraU(v,u*) e V@V*. (6) 

Then let cti, a n be a basis of V and a\, a* be a basis of V* such that (a* \ctj) = 
hj = 1) ■••) n ) where 5 is the Kronecker symbol. As a basis of V © V* one can 
choose (a\, 0)..., 0), (if) a*, a*), (tp a*, a*), then the matrix associated to the 
bilinear form in ([5]) becomes 

On In 
In On 

where n is the n x n zero matrix and I n is the n x n identity matrix. Accordingly, 
the basis given by 

Hi = ^ ((ai,0) + (ipa*i,a*)) , 
Xi = ((ai,0) - (^a*,a*)) 
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gives the diagonal form 

' In 



"In 

Then it is concluded that the bilinear form in ([5]) has signature {+1, +1, — 1, — 1}. 

n times n times 

□ 

Proof of Theorem [7J Existence of the map £ goes along with the lines with the usual 
argument [1] as follows. For every point x = (v,v*) £ V © V* the linear map £ is 
defined by £v* = v\ p *(p)-> where v is considered as a linear map on V*. To show that 
it is well-defined, suppose we have x = {v\,v\),x' = (^2^2) ^ ^ with v\ = v%- Then 
clearly (v± — V2,0) £ V which implies that v\ — V2 6 T> fl V. By the condition ([3]), 
this is equivalent to saying that [v\ — V2)\ p *(p) = ^ oi £v\ = £v^, as desired. The 
claim that the kernel of £ being T> (~)V* is trivial. 

Conversely, suppose that a subspace E C V and a linear map £ : E° — >• V/E are 
given. Then, if V is defined as in ((3D, it is easily seen that p*(T>) C (VnV)° is satisfied 
automatically. It remains to show T> has dimension n to satisfy the conditions of (ii) 
of Proposition (21 Since V includes subspaces E © {0} and C(E°) © E°, it is of at 
least n-dimensional. Observe that C can be extended to whole V*, and this gives an 
symmetric bilinear form on V* , defined by 

^{vl^):= l -{(vl\Cvl) + (vl\Cvl)). (7) 

Then by Lemma (H T> is a maximally isotropic subspace for an inner product on 
V © V* of split sign. Evantually dimD = n. □ 

Remark 1. If the linear map C in associated to a LD structure T> is skew- symmetric, 
then the LD structure turns out to be a Dirac structure. But in general C can be 
uniquely written as a sum of a skew- symmetric form L~ and symmetric form C + : 

£ = £- + £+. (8) 

If C is purely symmetric or skew- symmetric, then additionally one has 

P (v)° = vnv*. (9) 

To see this one can check that if £ is skew-symmetric then 

( P *(v) | pip n v)) = — ( P *(v n V) | pip)) = o, 

if £ is symmetric then 

(p*(D) | P (v n V)) = (p*p n V) | P (v)) = o. 

So, pip) C ip n V)° in both cases. A dimension count gives the equality, since 
dimip) = n. 

The equations ^ and are called the characteristic equations JIJ/. But the 
converse of the above claim is not generally true, that is, (0j and (0|) are not sufficient 
for £ to be symmetric or skew- symmetric. For instance, if £ is defined on whole V* 
and it is of full rank, then the characteristic equations are satisfied. 
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From the proof of Theorem [TJ LD manifolds have the following property. 

Corollary 1. If D is a LD structure on a vector space V , then there exists an inner 
product of split sign on V © V* which respect to which T> is a Lagrangian subspace. 

Theorem Q] gives a representation of a LD structure, another representation can 
be given as follows pQ. 

Theorem 2. Let T> be a LD structure on an n- dimensional vector space V , then 
there exist two linear maps a : M. n — > V and b : R n — > V* such that 

keradkerb = {0}, (10) 
(a(/cerb))° C rangeh. (11) 

Conversely, any LD structure T>onV can be defined by 

V = {(ax,bx); x G R n } C V @V* (12) 

for two linear maps a : 1" — > V and b : M n — > V* satisfying il(J\) and 

Proof. Let P be a LD structure on V. If one chooses a basis for D, then this is 
equivalent to giving two linear maps a : M n — > V and b : W 1 — > V* such that the 
basis becomes (ae 1; b ei), (ae n , b e n ). Since V is n-dimensional, (fTUj) is satisfied. 
Observe that T> n V = ker b and p*(T>) = range b. Then by the defining property ([3]) 
of T>, the relation (jlip is satisfied. 

Conversely, assume that P is given by (|12p . then by (|10p it is n-dimensional, and 
(fTT|) implies ([3]) by Proposition [2] (ii). □ 

2.2 Leibniz-Dirac structures on manifolds 

Let M be a n-dimensional smooth manifold. Consider a subbundle T> of the Pontrya- 
gin bundle TM © T*M. The following associated distributions and co-distributions 
to T> will be of importance throughout: 

G = {X £ TM; (X,0) G 2?}, 

d = {X G TM; (X,a) G 2? for some a G T*M}, 

and 

P = {a G T*M; (a, 0) G £>}, 

= { a G T*M; (X, a) G V for some X G TM}. 

For simplicity we will assume Go, G\, Po and Pi to be of constant rank. Recall that 
points on the distributions Go, G\ and the co-distributions Po, P\ have constant rank 
are called regular points. The set of regular points is an open dense set [TJ. 

If there are more then one LD structures, we will denote the bundles associated 
to V by Gq(V), etc. For any distribution G the co-distribution ann(G) is defined by 



ann(G) = {a G T*M ; (a\X)=0 for all X G G}, 
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and for any co-distribution P the distribution ker(P) is defined by 

ker(P) = {X £ TM; {a | X) = for all a £ P}, 

where ( | ) denotes the natural pairing between 1-forms and vector fields. 
Definition Q] can be given on a manifold as the following. 

Definition 2. Let T> be a smooth vector subbundle ofTM® T*M . Then T> is called 
a Leibniz-Dirac structure (LD structure) if 

Pi = ami G . (13) 

The following result gives a representation of a smooth LD structure whose proof 
is analogues to the linear case when considered pointwise. 

Theorem 3. IfDisa LD structure on an n- dimensional manifold, then there exists 
a vector bundle map C : Pi — > TM / Gq for which 

V = {(X, a) : X = Ca (mod G )}. (14) 

Conversely, if a distribution Gq on M of constraint rank, and a vector bundle map 
C : annGo — > TM/Gq are given, then a LD structure can be defined by |i^| ). 

Example 1. If the bundle map C is skew- symmetric, then the LD structure turns out 
to be a Dirac structure. This case includes many geometries like symplectic, Poisson, 
foliation, etc. 

Example 2. // there is a pseudo-Riemannian metric g on M , then it gives a bundle 
map Q : TM — > T*M. Then the inverse, say, C of this map gives rise to a LD struc- 
ture which is symmetric. This setting allows one to define gradient control systems 
J?}'. Note that one can also define a pseudo-Riemannian manifold as a Lagrangian 
subbundle with respect to the skew-symmetric pairing 

((X 1 ,ai) , (X 2 ,a 2 ))- = (ai | X 2 ) - (a 2 | X x ) 

for all (X±, ai), (X 2 , a 2 ) G TM © T*M which was given originally in JI]/. This will 
be studied in a future paper F73]/ . 

Example 3. In the general case, if C is as in Theorem^ and Pi = T*M, then LD 
manifold is simply called Leibniz manifold. These manifolds are studied in f^j but 
our definition also allows to add some constraints when modeling physical systems by 
LD structures (cf. Section^). 

As a linear LD structure, a smooth LD structure has the following property. 

Corollary 2. // a manifold M has a LD structure, then D is a Lagrangian subbundle 
with respect to a symmetric pairing on TM ®T*M . 
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Remark 2. As in the linear case C can be uniquely written as a sum of a skew- 
symmetric form C~ and symmetric form C + : 

C = C-+C + . (15) 

Observe also that C can be extended locally to a bundle map from T*M into TM which 
we also denote by C. Thus one obtains associated tensor fields B L : T*M xT*M — > R 
and B R : T*M x T*M ^ R by 

B L (ai,a 2 ) = (a 2 \Cai) (16) 
B R ( ai ,a 2 ) = (a 1 \Ca 2 ) (17) 

for all ai, a 2 E T*M. B L and B R are called left and right Leibniz tensors J2J/, respec- 
tively. Recall that there corresponds Leibniz parenthesis [, ] L and [, ] R associated to 
both Leibniz structures B L and B R , respectively, given by 

[f,g] L = B L (df,dg) and [f,g] R = B R (df,dg). (18) 
In the next section we will define an associated bracket to these tensors. 

Integrability of LD structures on manifolds is defined in accordance with the one 
on Dirac structures [T]. 

Definition 3. Let D be a LD structure on a manifold M . If the characteristic 
distribution Gq is integrable, then T> is called an integrable LD structure. 

In fact, this is a conclusion of a Dirac structure being integrable but in this 
generalized case of LD structures, integrability definition makes sense. This can be 
seen by the following result. 

Proposition 3. LetD be an integrable LD structure on a manifold M . If the foliation 
of Go is denoted by then M/& inherits two Leibniz structures, namely, left and 
right Leibniz structures. 

Proof. 1-forms on M/<3? are considered as ^-invariant 1-forms a € O x (M): a(T$) = 
0. Then induced structures are defined by means of ^-invariant 1-forms. They are 
apparently linear and well-defined, and so give Leibniz structures. □ 

Remark 3. Recall that, in case of Dirac structures the induced Leibniz structures 
coincide and they are specifically Poisson structures 

3 Dynamics on LD manifolds 

Dynamic properties of LD structures are given in this section. We first give the 
notion of admissible functions. The main ingredient of this section is a formulation 
of dissipative Hamiltonian systems with constraints. 
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3.1 Admissible functions 

Admissible functions on LD manifolds are denned as in the Dirac case [TJ. Let T> be 
a LD structure on a manifold M, and let T(T>) denote the smooths sections of T>. A 
function on M is called an admissible function if df G T(T>). If / is an admissible 
function, then (Xf,df) G T> for some vector field Xf on M. Accordingly, two brackets 
on admissible functions on M can be defined by 

[f,g] L =X f (g) (19) 

and 

[f,g] R = X g (f). (20) 

As usual let £ denote a representation of V and let B L and B R denote the left and 
right LD tensors respectively (cf. Remark [2]). Then 

[f,g] L = B L (X f ,X g ) (21) 

and 

[f,g] R = B R (X f ,X g ). (22) 

The following result is an extension of the Dirac case [I]. We only consider left 
Leibniz bracket, the same result follows for the right Leibniz bracket in a similar way. 

Proposition 4. With the notation above, admissible functions satisfy the Leibniz 
identity with respect to the left Leibniz bracket: 

[fg,h] L = g[f,h] L + f[g,h] L . (23) 

Furthermore, if f and g are admissible functions, so is the product fg. 

Proof. Let (Xf,df) and (X g ,dg) be two sections of P. Then 

[fg,h] L = X fg (h)=gX f (h)+fX g (h) 
= g[f,h] L + f[g,h] L 

which proves the Leibniz identity. To prove the other claim one computes 

(X fg ,d(fg)) = (gX f + fX g ,gdf + fdg) 

= g(X f ,df) + f(x g ,dg)ET(V). 

So, fg is an admissible function. □ 

From the above result, admissible functions on M form a Leibniz bracket which 
was defined in [9] not to be confused with Leibniz algebra brackets |16j . 

Remark 4. Recall the fact that the bracket of two admissible functions is again an 
admissible function, which follows from the Jacobi identity. Nevertheless, this may 
not hold in a manifold with a general LD structure. 



Remark 5. Observe that, a proof of Proposition^ can be given in terms of admissible 
functions as originally given in the Dirac case 
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3.2 Dissipative Hamiltonian systems with constraints 

Let T> be a LD structure on a manifold M, then one can extend the notion of implicit 
Hamiltonian systems [17J to (M,T>) as follows. 

Definition 4. Let H : M — >• R 6e a Hamiltonian. The dissipative implicit Hamilto- 
nian system (DIHS) corresponding to (M,D,H) is given by 

(x,dH(x)) G £>(x), x G M. (24) 

In this setting, Gi describes the set of admissible flows and Pi describes the set 
of algebraic constraints. 

If the splitting (fT5|) is considered, by ([26]) one obtains 

This equation has several meanings by the following specific examples. 

Example 4. If the bundle map C is skew- symmetric, then the GIHS in \2J$ is an 
implicit Hamiltonian system and the equation H25]) gives the conservation of energy. 
This special case of conservative Hamiltonian systems with constraints have been 
studied by many authors in Dirac structures framework 1181 11 , 0, {(Jjj. 



Example 5. If C is symmetric then the system \2J$ is a generalization of pseudo- 
gradient systems [Tjj. This can be seen as follows. If £ is taken as the restriction 
to TM/Gq of inverse of the musical isomorphism g* : TM — > T*M of a pseudo- 
Riemannian metric g on M , then one obtains a gradient system with constraints. In 
that case the equation 1125]) is called the entropy equation /?/. As was pointed out in 
Exapmle it is more advantageous to handle these systems in terms of Beltrami- 
Dirac structures f!5f . 

Example 6. If C is arbitrary, that is, not necessarily symmetric or skew-symmetric, 
then the GIHS is a generalization of metriplectic systems with constraints \19§ . In 
this case the equation 125]) is may be considered as the dissipation equation [19]. This 
general case is used to model so many physical systems. For example, systems with 
nonholonomic constraints, dissipation, etc JEH EH [23 EE EH EH EOF- 

Remark 6. Assume that T> is represented in the form |i^[ ), then the DIHS corre- 
sponding to (M,V,H) has a local representation 

x = C(x) ^-(x) + g(x) A, 
= g T (x)^-(x), 

where ^r{x) stands for the column vector of partial derivatives of H , and g(x) is a 
full rank matrix with Img(x) = G$(x), and A are Lagrange multipliers corresponding 
to the algebraic constraints = g T (x) ^jr(sc) \17$ - 
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Example 7. Let Q be an n- dimensional configuration space and let q = (qi, 
be a local coordinate system on Q. Kinematic constraints on Q are given by 

A T (q)q = 0, 



(27) 



where A(q) is an n x k matrix, k < n, whose rank is k for all q E Q. Then the 
constrained dissipative Hamiltonian system on T*Q is defined by equations 



q 




o n 


. p _ 




. "In 



In 

R(g) 







+ 


O n 






. Mq) . 



A. 



Hfe P) 



1 

dp 



9h -(q,p) 



= [ O n A T (q) 
where R(g) is a semidefinite matrix. 

Example 8. (\2$ ) We consider a particle moving on a surface z = yx in M 3 with a 
friction force proportional to the particle velocity. The Hamiltonian is given in terms 
of cartesian coordinates x,y,z and their conjugate momenta by 



H(x,y,z,p x ,p y .p z 



(28) 



and the characteristic distributions are given by 
G = span{^- -y-^—} 

OPz OPx 

P\ = span{dx,dy,dz,ydp z + dp x ,dp y }. 
The corresponding LD structure C in matrix form is given by 

C = 3 - R, 



where 



R 



o 3 

-Is 

o 3 



I 3 

o 3 



Ml 






o 3 











with jUj(g) > is the directional and space- dependent damping coefficient fWjj. Then 
the equations of motion read 
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4 Reduction and transformation of LD struc- 
tures 

In this section, we first discuss reduction of LD structures on subspaces and more 
generally on subbundles. Then in Subsection 14.2} we consider maps which push 
forward LD structures. 

4.1 Reduced LD structures 

Let D be a LD structure on a vector space V. If £ C V* is a subspace, then it is 
possible to reduce the LD structure on S. For this purpose, it is desirable to use the 
representation of the LD structure given in Theorem [TJ 

If T> is given by a subspace E C V a linear map £ : E° —> V/E, then the reduced 
LD structure Ps is defined by the subspace E% = ker(E° n S) and the restriction 
£ s of £ to E° n S. Then £> s is given by 

p s = {(w,w*); w = £ E w*(mod£ s )} CS*®S. (29) 

Observe that an alternative way of defining reduced Dirac structures on subspaces 
of V is not trivial in general. Because, this makes sense if linear maps of form 
oj : V —> V* exist. However, if £ is symmetric or skew-symmetric one can also define 
a linear map from a proper subset of V into V* and the procedure of reduction given 
in [1] works. 

The construction of reduction given in the linear case can also be considered for 
smooth LD structures. If P is a smooth LD structure on a manifold M, and if 
£ C T*M is a smooth submanifold, then the following clean intersection property 
is required to reduce V. Let us denote the reduced LD structure by V%. As the 
characteristic distribution of £>£ will be P\(T>)r\T* M , this distribution should 

be a subbundle. Then one can define the reduced LD structure as an extension of 

422). 

4.2 LD maps 

In this subsection, we study maps which carry LD structures to LD structures. We 
start with the linear case then consider smooth LD maps. Let Dbea LD structure on 
a vector space V. As the LD structure V is associated to a linear map £ : V* — > V, 
we can be pushed forward T> into another vector space by a linear map |14j . 

Definition 5. Let eft : V — > W be a linear map between vector spaces, and let T>y be 
a LD structure on V . Then push forward image ofDy is defined by 

T^Vy) = {(0v,w*); v G V,w* e W*,(v,</>* w*) G V v } dW ®W\ (30) 

where <ft* : W* —> V* denotes the adjoint map of <ft. 

Proposition 5. Let (j) : V —> W be a linear map. If V is a LD structure on V , then 
TcpiDy) is a LD structure on W . 
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Proof. Let us denote the subspace F^Vy) by T>w-> an d let pv,p v (resp. pWiPw) 
denote the projections from V ® V* (resp. W © W*) onto the first and the second 
factors respectively. First notice that 

V w nW = (f){V v C\V). (31) 

Let us consider the associated linear map Cy : P\(Vy) C V* — > V to Vy which is 
defined in Theorem [TJ Then there is a well-defined map Cy/ from p^(F ( f > (Vy)) into 
W given by 

C w w* := cj>(C v ((j>* w*)). (32) 
We claim that F<f>(Vy) = Vyr where 

V w = {(w,w*); w = C w w* (mod D w nW)}. (33) 

which is then a LD structure on W by Theorem [TJ We first show that F^Vy) C Vy/. 
In fact, (w,w*) G F^Vy) implies w = 4>v with (v,cj)* w*) G Vy, then 

v = C y {(f)* w* ) (mod (V D V)) (34) 

by the definition. This shows, by considering the images of both sides of the equation 
(j3U) under 0, that 

u> = £vv ^* (mod D\y n VF) 

by (f3"T|) and ([32]) . Then (w,w*) G D^. To show that £V C F^iVy) we assume that 
(w,w*) G XV > namely w = £ w u>*(mod£V n W). This implies by ([3TJ) and (f32|) 
that w = c/>(£y(</>* ^*)) + for some w' G XV H V. Since is linear u; = </>t> for 
v = £y(</>* «>*) + w' which implies that (v, 4>* w*) G Vy. Then J^(Py) C XV- This 
concludes the proof. □ 

By the proof of the preceding result the following conclusion is obtained as in the 
Dirac case [T4"] . 

Corollary 3. Let <fi : V —> W be a linear map, and T> be a LD structure on V , Then 
for the Dirac structure Vy/ = F^iVy) on W the following equalities hold: 

v w nw = <p{Vy n V) (35) 

and 

Cy/ = <fi o Cy o (j)* . (36) 

Remark 7. Alternatively, one can follow the methods given in to study maps 
between vector spaces with LD structures. As indicated in the setting for the 
Lagrangian relations ]25^ is valid between inner products of split sign onV ®V* and 
W © W* . By using Corollary one can set suitable inner products on both spaces 
and an alternative proof of Proposition [5| can be given. 
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Now we proceed with the smooth case. Let M be a manifold and let T>m C 
TM © T*M be a LD structure. If <p '■ M — > N is a smooth map then the forward 
image of T>m is defined by 

^(%) = {(MX),P); (X,<f>*{p)) G (37) 

Clearly, J^Pm) may not be a subbundle, to ensure its being a subbundle the fol- 
lowing definition is given. Dm is called ^-invariant if 

T^Pm)* = {(M*),P) ■ (X><t>*(P)) G (^m)*} (38) 

is invariant for all x on the same <^-fibre [26]. Nevertheless, <^>-invariance is not 
sufficient, one also needs some extra assumptions as in the following [26J. 

Proposition 6. Let <f> : M —> N be a surjective submersion and Dm be a LD 
structure on M. If ker {(/)*) V\Dm has constant rank and Dm is ^-invariant, then the 
forward image ofT>M defines a LD structure on N. 

Proof. First, denote T^,(Dm) by Dn, then Gq(Dn) = 4>*(Go(T> m )). As a result the 
rank of Gq(T>n) is constant if and only if the rank of ker(<^*) n Dm is constant, since 
Gq(Dm) is a subbundle by the hypothesis. Thus Dn is a subbundle. Next, by the 
0-invariance, one can define a tensor field C : T*iV — > TN by 

C N (<t>(z)) v ■= <j>{C M {xWr,)). (39) 

Then J-^Dm) can be defined in terms of Cn as in Proposition □ 

Remark 8. As it is pointed out in }26j . if there is a Lie group action on M , which 
preserves the LD structure, then the conditions in the hypothesis of Proposition [7J 
are satisfied if furthermore the action is free and proper. This and related questions 
about reduction will be the subject of a future work. 

5 Conclusions 

We defined linear and smooth Leibniz-Dirac structures which are generalizations of 
Dirac structures. We studied the geometry and dynamics of LD structures both 
on linear spaces and manifolds. It was explained with several examples that LD 
structures are capable of formulating dissipative implicit Hamiltonian systems with 
constraints. We hope that LD structures will find applications in physics and related 
area. 

However, there remain many questions on both geometric and dynamics proper- 
ties of LD structures unexplored, which we address herewith. As it is known, a more 
general setting of Dirac structures on Courant algebroids has a growing importance 
\27\ 128] . Accordingly, LD structures may be defined on generalizations of Courant 
algebroids such as Leinbniz algebroids [22 [3U]. As also pointed out in Remark [SI 
one can study reduction of symmetric LD structures as a generalization of Dirac 
reduction [281 13 lj . Another related problem is quantization of dissipative systems 
[32j . We hope our results will contribute to this and related problems on dissipative 
Hamiltonian systems with constraints. 
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